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snsoDDcncN 


A technique for electronically rotating the polarization basis of a 
orthogonal-! i near polarization radiometer has been successfully 
demonstrated by investigators at the Georgia Institute of Technology. Ihe 
technique is based on the measurement of the first three feedhorn Stokes 
parameters along with subsequent transformation of this measured Stokes 
vector into a rotated coordinate frame (Gasiewski, 1990; Parks, et al, 
1988; see also Appendix A) . The technique requires an accurate measurement 
of the cross-correlation between the two orthogonal feedhorn modes, for 
which an innovative polarized calibration load has been developed . 

Ihe electronic technique can simplify designs of satellite-based 
scanning polarized microwave imagers by electronically (rather than 
mechanically) rotating the polarization basis of the radiometer. For 
scanning imagers , some means of rotation is required in order that the 
feedhorn basis remain coincident with the Earth's vertical-horizontal basis 
throughout the entire scan. Potential applications include the Defense 
Meteorological Satellite Program (DMSP) microwave precipitation 
instruments, the TFMM Microwave Imager (TMI) , the EOS Multi frequency 
Imaging Microwave Radiometer (MIMR) , and airborne passive microwave 
radiometers for meteorological research. 

The experimented portion of this investigation has consisted of a 
proof -of-concept demonstration of the technique of electronic polarization- 
basis rotation (EPBR) using a ground-based 90-GHz dual orthogonal-linear 
polarization radiometer. I>ractical calibration algorithms for ground-, 
aircraft-, and space-based instruments have been identified and tested. 

Ihe theoretical effort has consisted of radiative transfer modeling using 
the planar— stratified numerical model desc r i bed in Gasiewski and S tael in 
(1990) . 
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SUttEFY OF ACTIVITIES 


Activities sponsored by this grant at Georgia Tech during the period 
from June 1, 1990 to November 30, 1991 include acquisition of required 
hardware for the ground-based EFBR demonstration instrument, construction 
of all electri cal and mechanical subsystems for the ground-based 
demonstration unit, development of practiced calibration algorithms, 
ground-based experiments verifying the concept of EFBR, and radiative 
transfer calculations. 

a. Hardware 

Briefly, the demonstration instrument consists of a 90-GHz dual 
orthogonal-linear mode radiometer, modified to include an additional cross- 
correlating channel (Figure 1) . The addition of this channel allows the 
orthogonal-linear polarization basis of the radiometer to be rotated by any 
desired angle. Hie radiometer is supplied by the NASA Goddard Space Flight 
Center, and will be available for use at Georgia Tech through March, 1992, 
and beginning again in June, 1992. 

The cross-correlation channel is implemented as an adding correlator. 
Relatively little hardware beyond that required to implement a conventional 
dual-polarization radiometer (i.e. , without cross-correlation channel) is 
required. This hardware consists of three microwave power 
divider/ccmbiners, a video detector and a video summing amplifier. In order 
to meas ure the cross-correlating signal with minimal corruption by additive 
radiometric noise, the total gains (including IF, RF, and video) of all 
channels must be balanced, and the radiometric noise from the vertical and 
horizontal channels mus t be subtracted in the sunming amplifier (Gasiewski 
1990) . 

An adjustable local oscillator (W) phase shifter was installed in one 
of the two LO-mixer paths to allow precise adjustment of the cross- 
correlation channel to measure either in-phase correlation (Ty) or 
quadrature correlation (Ty) . The adjustable phase shifter provides the 
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capability to meas ure all four Sto k es parameters by switching the ID phase 
by 90°. 

An IEM PS-2 Model 30 computer with a Keithley Metrabyte A/D converter, 
digital I/O system, and stepper motor drive capability was configured to 
provide temperature and radiometric data acquisition, instrument control, 
and calibration capabilities. Hie data acquisition and calibration software 
were developed at Georgia Tech. 

The measured responses of the three channels to a li ne a r ly polarized CVJ 
source as a function of source polarization angle are shewn in Figure 2. 
The data indicate good cross-correlation response, as seen by the 45-degree 
phase shift of the U-channel response relative to the A- and B-channels. 
Worst-case cross-polarization isolation for the A- and B-channels is 
approximately 15 dB. Although this is somewhat low, it appears to be 
acceptable for EPBR. (It is noted that a 15 degree offset in the electric 
field orientation relative to the instrument's window mounting plate at the 
co— polarized response maximum of the A- channel was observed. This is 
suspected to be caused by imperfect alignment of the radiometer's 
orthogonal-mode coupler. The misalignment is not serious for this 
investigation, since it remains constant. However, it should be accounted 
for by other investigators using the GSFC 90-GHz instrument.) 

B. Calibration Algorithms 

The radiometric response of all three channels must be known in order 
to interpret the measured output voltages. A polarized calibration load 
using hot and cold microwave absorbers and a polarization-splitting wire 
grid (Gasiewski, 1990; Figure 3) was constructed for this purpose. The 
configuration of the absorbers and grid were chosen to provide a precisely- 
known polarized Stokes field to the feedhom of the radiometer. 
Observations of this load at several angles of rotation allowed the gains 
and offsets of all three channels (both orthogonal mode and cross- 
correlation) to be measured. 

In calibrating the radiometer, four parameters must be determined for 
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each of the three channels: a primary gain, an offset , and two crosstalk 
gains. Althou^i the crosstalk gains should be zero, in practice they are 
not and need to be determined for accurate EPBR. To determine the four 
unknown parameters, an invertible gain/ offset observation matrix must be 
obtained for each channel. For example, let the output voltages of the U- 
channel for N calibration looks be given by the values v^, and let the 
corresponding feedhom Stokes' brightnesses be (T^ir Tgi/ ^Ui) • These 
values are linearly related by the U-channel gains 9uB' 9uu) arK ^ 

offset Oy: 


where nj_ is a radiometric noise voltage. If the above matrix is rank four 
(i.e. , full row rank) and the noise standard deviation is small, then a 
unique solution for the gains and offset exists. In principle, the 
radiometric noise can be reduced arbitrarily by increasing the integration 
tim e of the calibration look or (equivalently) averaging over a number of 
looks. 

A rank-four observation matrix can be obtained using four calibration 
looks: three at the polarized load rotated to three appropriate angles 
(e.g., 0,45, and 90°) and one at an additional unpolarized load. This 
scheme was found to most convenient for laboratory experiments, and was 
adopted. The additional view of an unpolarized load was accomplished by 
sliding a piece of rocm-temperature absorbing foam over the feedhom. 


The expected accuracies of several alternate schemes for calibrating 
the instrument were also investigated. The performance of each of these 
schemes was analyzed using the matrix software package MATIAB. Each scheme 
consists of a particular hardware configuration employing a polarized 
calibration load and/or an unpolarized load viewed in a particular 
calibration load viewing sequence. A calibration sequence consists of a 
series of views of either an unpolarized load or a polarized load rotated 
to predetermined angle. 
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Indeed, the best simulated results cure achieved using a calibration 
scheme with a rank- four gain/ offset observation matr ix. Non full - rank 
observation matrices and/or noise-corrupted calibration looks are less 
desirable for measuring the channel gains and offsets. For exairple, if the 
calibration sequence is implemented by viewing only a simple two- 
temperature polarized load, the crosstalk gains can not be unambiguously 
determined. This is not to say that the crosstalk gains could not be 
statistically determined, but this would require a-priori statistical 
knowledge of the drift behavior of the hardware. Thus, the non-full rank 
schemes are highly specific to the radianeter hardware, and not generally 
recommended. 


Thus, for airborne or spacebome applications, a view of an unpolarized 
load (e.g. , cold space) along with the simple two-temperature polarized 
load during the calibration cycle is nearly essential. However, a slightly 
more complicated three-temperature polarized load can provide nearly the 
same results without the necessity of viewing an unpolarized load. Only a 
polarizing wire grid need be ro tated , thus minimizing mechanical 
complexity. (A splash plate would also have to be inserted to allow the 
feedhom to view the calibration load.) This scheme would use one hot 
calibration load (viewed by the polarization transmitted throu^i the wire 
grid) and two cold calibration loads (each viewed by the polarization 
reflected from the wire grid) . One of these two cold loads would typically 
be an ambient-tenperature blackbody (approximately 250 K) , and the other 
could be cold space (2.73 K) . A calibration sequence viewing both the 
hot/ambient and hot/space load combinations at appropriate angles of the 
rotating wire grid provide a full-rank gain/offset observation matrix. 
Since only the wire grid need be rotated, the rotating mass is minimal, and 
slip rings for electrical connections are unnecessary. 

C. Ground-based Demonstra tion of EPBR 

The ability to calibrate an EPBR radiometer and transform the feedhom 
brightness temperatures into the natural basis was successfully 
demonstrated using outdoor observations of a reflecting water surface at 
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near-Brewster angle (65°) incidence under clear-sky conditions . For the 
chosen incidence angle , the vertically and horizontally polarized 
brightnesses differ significantly (Ty 260 K, T^j 135 K) , causing a high 
degree of polarization (m'0.32) . A conically-scanning mirror was configured 
to sweep the radiometer beam across the water surface, thereby rotating the 
feedhom polarization basis relative to the natural basis while maintaining 
constant angle of incidence. 

Results of a water surface scan are illustrated in Figure 4, along with 
predicted brightn ess es based on numerical radiative transfer calculations 
(see also Appendix B) . The transformed feedhom brightnesses track the 
predicted values well, as indicated by the nearly-constant values of Ty, 
Tr, and U. Ideally, these would remain constant throughout the scan. 
However, a small systematic error of less them 7 K was caused by mechanical 
misalignment of the scan mirror and feedhom axes. In practice, precise 
alignment would eliminate this error. Other sources of discrepancy between 
the predicted and measured brightnesses include errors in the assumed 
atmospheric state, which was taken to be the US s tandar d atmosphere for the 
purpose of this comparison. 

D. Alternatives 

As an alternative to EPBR and the associated calibration scheme, 
another means of achieving electronic rotation of two orthogonal-linear 
feedhom modes was considered. This technigue uses ferromagnetic devices 
based on the Faraday effect. These devices can provide selectable rotation 
of orthog 

onally-polarized circular waveguide modes in response to the strength of an 
applied DC magnetic field. However, they have the following shortcomings 
for radiometric applications: 

* The bandwidths are inherently small (1%) • At least 2 GHz 
bandwidths are desirable for radiometric observations near 
90 GHz. 

* They are somewhat lossy (>1 dB at 90 GHz) , and hence reduce 
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radiometer sensitivity. 


* Practical Faraday rotators are currently unavailable for 
frequencies above 100 GHz. (Future passive microwave observations of 
cirrus observations may require polarized observations at frequencies 
as high as 340 GHz.) 

* At best, Faraday rotators can be used only over a single radiometric 
bard. This precludes their use in multibard feed systems, such as 
SSM/I or TMI. 

The EFBR technique overcomes all these limitations with minimal 
hardware complexity and the added capability to observe both the third and 
fourth Stokes parameters U and V. 

E. Radiative Transfer Model 

Most of our efforts during the period of this grant have focussed on 
the construction of hardware for the ground— based demonstration instrument. 
However, seme progress was made in upgrading the Georgia Tech MPT numerical 
radiative transfer model. This code, the original version of which was 
written by the author at the Massachusetts Institute of Technology in 1987 , 
is a comprehensive scattering— based numerical radiative transfer model. 
Currently, the model incorporates scattering and absorption from liquid and 
frozen spherical hydrometeors. In order to interpret polarized microwave 
brightness temperatures observed over precipitation, the model is being 
upgraded to accommodate aspherical liquid and ice particles. 

The MPT numerical model uses an iterative (or perturbation) approach to 
calculate the temperature weighting functions for a planar-stratified 
scattering and absorbing atmosphere. Hydrcmeteors are currently modelled by 
Marshall-Palmer distributed liquid spheres and Sekhon-Srivastava 
distributed ice spheres. To be useful in interpreting polarized microwave 
brightness temperatures caused by aspherical hydrcmeteors, modifications to 
the phase matrix, extinction matrix and absorption vector computation are 
required. Development of subroutines to calculate these quantities for 
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planar-stratif ied distributions of prolate and oblate spheroidal 
hydrcmeteors is ongoing. 

For electrically small particles, the Raylei^i-Gans approximation will 
be used. For electrically large particles, an exact solution using 
spheroidal coordinates for the scattered field of a single spheroidal 
particle (Asano and Sato, 1980) was investigated. However, two attributes 
make this approach undesirable: (1) lack of reliable routines for 
calculating the vector spheroidal wave functions required by this 
formulation, and (2) the number of nested sumnations required. The T-matrix 
approach (Tsang, et al, 1985) is more desirable because the formulation is 
based on the more familiar vector spherical modes, and the number of nested 
summations will be fewer. 

A major problem in aspherical hydrometeor radiative transfer modelling 
is determining realistic models of ice hydrometeor particle size, shape, 
and orientation distributions, particularly within cirrus anvils and 
thunderstorm tops . To this end, we are investigating a statistical 
relationship between aspherical particle orientation and the local 
atmospheric electric field. In order to minimize electrostatic energy, 
particles will preferentially align themselves with their longest axis in 
the direction of the field (Stratton, 1941). The resulting canting angle 
distribution can be ccnputed using a Boltzman probability distribution. 


The rel ationship between the local field and particle orientation has 
been hypothesized to be the cause of observed changes in the phase delay 
and depolarization of microwave signals received from satellite beacons 
(Cox and Arnold, 1979) . It is also plausible that polarized microwave 
brightness temperatures could be used to identify regions of atmospheric 
electrification, or vice-versa. 
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Figure 2. Responses of auto- and cross-correlation channels to a 90-Giz CW 
source as a function of the polarization angle /p of the source. 
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Figure 4. Results of near-Brewster angle conical scans of a highly 
polarizing water surface. The parameters T^, Tjj/ and Tg are the measured 
feedhom brightnesses; Ty, %, and U are the resulting transformed 
brightnesses in the Earth's vertical-horizontal basis. The brightnesses 
predicted using a numerical radiative transfer model are illustrated by the 
heavy lines. 






CCNCLDSICTE AND REOCMfflNDftTICKS 


Hie results of the calibration experiments and water scans effectively 
demonstrate that EPBR is a practical alternative to mechanical rotation of 
a dual -polarization radiometer about its feedhom axis. Applications of the 
technique include both conical and cross-track scanning instruments , 
although the cross— track scanner benefits less since polarization 
information is typi cal ly not useful at observation angles less than 30° 
from nadir. For an instrument designed to measure all four Stokes 
parameters, the EPBR technique has a decided advantage, since switching 
between measurements of U and V requires only a single ID phase shift. 

Our recommendations with respect to EPBR and its applications are 
fourfold: 

(1) It would be interesting to ccmmence ground-based observations of 
polarized microwave emissions from precipitation cells. Such observations 
could be used to investigate correlations between oriented ice particles 
and thunderstorm electrification. To support this investigation, a two-axis 
imaging scanhead is being planned for use with the 90— GHz ground-based 
radiometer. To complement the passive measurements, coincident downlink 
scintillation and depolarization measurements using a satellite beacon 
(e.g. , ACTS) are being considered. The local electric field strength can be 
measured using an electric field mill probe. The results of coincident 
passive, active and electric field measurements could provide an improved 
understanding of propagation through regions of oriented ice particles, 
particularly thunderstorm anvils. 

(2) To provide polarized high resolution images of precipitation cells 
as would be observed from a satellite, the implementation of an aircraft- 
based imaging radiometer with EPBR is recommended. Two practical options 
might be considered. The first is a modification of the Advanced Microwave 
Precipitation Radiometer (AMFR) to include EPBR. (A proposal describing 
suggested EPBR modifications has recently been submitted by the authors of 
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this report to the contract monitor at NASJyMSPC. ) This would be a 
relatively inexpensive means of achieving broadband dual-polarization 
imaging capability, but suffers frcm two ncn-critical problems. First, the 
mechanics of retrofitting the AMFR with suitable polarized and unpolarized 
loads will be difficult due to the space restrictions. Second, the 
resulting cross-track scanned images would exhibit useful polarization 
information only near the ends of the scan. 

A second cption (albeit more expensive) is to construct a new 
conically-scanned instrument based on EPER for flights on either the NASA 
DC-8 or ER-2 . Suggested frequencies include window channels at 18, 37, 90, 
155, 220, and 340 GHz. The instrument wculd be particularly useful for 
studying polarized emissions frcm electrified thunderstorm anvils and thin 
clo uds over water backgrounds. By virtue of its conical-scan imaging mode, 
such an instrument could also provide accurate underflight data for the 
SSH/I, TMI, and MIMR instruments. 

(3) Meas urement and interpretation of the full Stokes vector over a 
wave-covered water surface at oblique viewing angles will provide 
information needed to optimize remote sensing techniques for ocean surface 
winds. Under appropriate conditions, it has been shewn to be possible to 
remotely sense the direction of propagating gravity or capillary waves on 
the water surface, thereby providing a measurement of surface wind 
direction (Wentz, 1991). 

To further investigate the potential abi l ity of a satellite-based 
conical ly— scanned radiemeter to remotely sense ocean surface wind 
direction, we recommend measurements using the 90-QIz ground-based 
instrument observing a simulated wave-covered ocean surface. Thus, the 
sensitivity of the Stokes' parameters Ty, T^j and V to wave direction can be 
investigated. All of these parameters can be measured simultaneously vising 
the 90-GHz ground-based radiometer provided that the ID phase at one of the 
two mixers is shifted by 90°. This can currently be done manually. 
Simulated ocean wave models of various heights and periods can be 
constructed from a lew— loss, low— dielectric constant matrix (s.g. , open- 
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cell f oam ) saturated and covered with oceanic s al ine. In this manner, the 
emissivity of an ocean wave in the geometric optics limit can be studied in 
detail as a function of observation angle. 

A means for electronically (as opposed to manually) shifting the ID 
phase in order to facilitate fast, reliable switching of a spacebome EPBR 
radiometer from the U measurement mode to the V mode should also be 
demonstrated. This can be accomplished using ferrite switched circulators. 

(4) In order to interpret aircraft- and ground-based thunderstorm and 
simulated ocean wave observations, it is recommended to continue the 
polarized radiative transfer modelling commenced during this grant. This 
modelling emphasis has two components: (A) millimeter-wave radiative 

transfer modelling within oriented needle- and plate-shaped ice crystals 
using an ellipsoidal particle model, and (B) reflectivity and emissivity 
modelling of lossy periodic surfaces, in particular, capillary and gravity 
ocean surface waves, viewed at oblique angles. 

Our specific recommendation is to upgrade the Georgia Tech MRT 
numerical code. The overall goal will be to use the model to interpret the 
experimental data, and to develop passive microwave retrieval algorithms 
for cloud and precipitation parameters, and ocean surface wind direction. 
Aspherical KT calculations will be essential in developing precipitation 
parameter retrieval methods using polarized microwave observations, and 
interpreting both airborne and ground-based polarized microwave 
observations. We also recommend investigation of the hypothesized 
relationship between polarized microwave brightness temperatures caused by 
frozen hydrcmeteors and the local atmospheric e l ectric field. An additional 
u se f u l capability would be full Stokes vector calculations (i.e. , all four 
parameters^ over a wave-covered ocean surface. This upgrade will be 
critical in investigating the ocean surface wind-direction sensing 
capabilities of full -polarization passive microwave brightness 
measurements. 
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Abstract- A combined hardw art and software technique is 
dweribed for the measurement of the vertically and horUontally 
polarised components of the atmoepbenc microwave brightness 
temperature at arbitrary scan angles using a fixed-feed dual- 
polarisation conical or croee track scanner. The t e chni que is an 
implementation of electronic polarization basis rotation, and is 
facilitated by the correlation of the in-phase amplitudes of two 
orthogonal linearly-polarised feedhorn modes. A linear minimum 
mean-squ are-error estimator is subsequently used to infer the ver- 
tical and horizontal brightnesses from the radiometer outputs. A 
practical hardware configuration and calibration scheme is pre- 
sented. 

The technique is particularly useful for passive microwave 
precipitation parameter retrieval using spacebornt or airborne 
conical scanners. The ultimate performance of the technique de- 
pends on the calibration of the receiver and balancing of the 
radiometer channels. Laboratory and aircraft flight testa of elec- 
tronic polarisation basis rotation are planned. 

Key Words - radiometer, polarisation, microwave, precipi- 
tation. 


1. Introduction 

Recent investigations into passive microwave remote sens- 
ing of precipitation have suggested that the difference between 
the vertically and horizontally polarised brightness temperatures 
contain useful information concerning the intensity of precipita- 
tion over ocean [1,2], and possibly, the presence of aspherical 
particles in the cell top [3]. The polarisation difference can also 
be used as an aid in detecting snow cover and discriminating 
between land and water backgrounds [4]. 

For the terrestrial troposphere and lower stratosphere, and 
for frequencies at which non-reciprocal propagation effects caused 
by the Earth’s magnetic field can be neglected, only the vertical 
and horisontal Stokes’ parameter brightness temperatures {Tv 
and Tjf) art non-sero. Hence, it is appropriate to design earth- 
probing instruments to measure, at most, these two parameters. 
Because the V - E polarisation basi# requires a miminal number 
of non-sero parameters (2) to uniquely describe the radiation 
field, it is termed the "natural” polarisation basis. 


2. Measurement of Tv and Ts 


Ideally, the polarisation basis of an airborne or spaetborne 
dual linearly-polarised radiometer should coincide with the natu- 
ral polarisation basis, so that the two radiometer channel outputs 
will be (upon proper calibration) Tv nud Tjj. Unfortunately, 
the need to image extended regions of the atmosphere by phys- 
ically scanning the antenna beam makes it difficult to build an 
instrument which retains polarisation coincidence with the nat- 
ural basis at every spot. Consider two common scanner config- 
urations for passive microwave imaging: conical and cross- track. 
The conical scanner is typically consists of an off-axis parabolic 
antenna rotating about the feedhorn axis, thus sweeping the an- 
tenna beam through a cone-shaped swath. The croee- track scan- 
ner consists of a rotating scan mirror oriented 45* with respect 
to the feedhorn axis. The resulting wedge- shaped swath is typi- 
cally oriented transverse to the motion of the platform, yielding 
a raster scan below the flight track. 


In their simplest configurations, both conical and cross- 
track scanners employ a feedhorn fixed rigidly to the instrument 
platform. For dual linearly-polarised feedborns, this results in 
the feedhorn polarisation basis rotating (with respect to the nat- 
ural basis) during the scan. The angular difference ^ between 
the two bases is termed polarization basis skew (the two bases 
coincide at 4 = 0). To obtain sero polarisation skew at all Kan 
angles, the feedhorn might be mechanically rotated along with 
the scanning reflector. However, alignment of the feedhorn and 
natural polarisation bases by mechanical rotation is both cum- 
bersome and expensive. 


As sufgested by Parks et si. [5] for the conical scanner, a 
more desireable way to eliminate the polarisation skew through- 
out the scan is to employ an inexpensive fixed-feed scanner, and 
to electronically rotate the polarisation basis. The linear trans- 
formation of brightness temperatures T from the natural polar- 
isation basis to feedhorn temperatures 7/ in the instrument's 
polarisation basis is: 


( Tf A \ eos*d «n*d Jsin2^ 

TfB I * | sin 2^ 

T/c J -sin 2^ sin 2# cos 2^ 


Tv 

Tb 

Tv 




0 ) 


where Tfu * R*(EaEb) *• the cross-correlation between the two 
orthogonal linearly polarised feedhorn modes (A and £), and d 
is related to the scan angle. In practice, the measured value of 
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7f (denoted ?/) is corrupted by additive instrument noise due 
to calbration errors and finite integration time: 

f f = V(+)-T + H r (2) 

With electronic polarisation basis rotation, the brightness 
temperatures in the natural basis are estimated using the mea- 
sured brightnesses in the feedhorn basis. The estimation of Ty 
and Tm are facilitated by the measurement of 7/p along with 
T/a and T /*. We addrees three issues concerning the practi- 
cal implementation of electronic polarisation basis rotation on 
both cross track and conical scanners: 1) the hardware config- 
uration of a 7/p -channel exhibiting acceptably small crosstalk 
and instrument noise. 2) calibration of the 7/p -channel, 3) an 
algorithm for estimating the brightness vector 7 in the natural 
basis from the measured data in the feedhorn basis. 

t Hhdaaa 

Fig. 1 illustrates a practical configuration of microwave 
and analog hardware for the measurement of the three compo- 
nents of 7/. The radiometer consists of a single superheterodyne 
down-conversion stage for each of two channels A and B, and 
a cross-correlating channel U. The phase shift is adjusted 
to cancel any RF path-length differences between the A and B 
channels. Square-law detectors and integrators provide voltage 
outputs proportional to the signal power in the respective inter- 
mediate frequency (IF) channels. 

lb implement the 27-channel, the A and B IF signals are 
added prior to detection. The output voltage vector V is thus 



* yT/ + o + u* (3) 

where J and Z are instrument gain and offset parameters, respec- 
tively, and fit, w J j* is the instrument noise referred to the video 
outputs. Upon calibration, the feedhorn brightness temperatures 
are subsequently found from the output voltages: 

ff * F’Mv-S) (4) 

Is order to detect the relatively small cross-correlation 
signal Tfjj in the presence of the larger signals 7/* and 7/$, the 
video outputs of the A and B channels are subtracted from the U- 
channel video output. Provided that the total IF and video gains 
associated with the power splitters, equare-law detectors, and 
integrators are matched to within — 0.5 dB, the gain parameters 
$av and fBU will be small compared to guv (thus minimizing 
crosstalk), and the instrument noise ft? will be nearly white, 
with covariance: 

^KrKr * m wf7 (5) 

where 7 is the identity matrix, and (•)* is the transpose operator. 
The noise standard deviation is: 

- - a* (,) 


where Tgri » the eyetem now* temperature (receiver plui in- 
tense), W is the IF bandwidth, end r ie the effective integration 
time. In practice, balancing is accomplished by video and IF 
component matching and thermal sta bilis ation. 

L Calibration 

Calibration of all channels can be accomplished using a 
polarised calibration device (Fig. 2) consisting of two highly- 
absorbing microwave black body loads and a polarisation-splitting 
wire grid. One load is cooled to a low kinetic temperature Tcold . 
and the other heated to a high temperature T*or« where T B ot - 
Tcold » 80 K. The wire grid combine* the two load bright- 
usm temperatures into two orthogonal- linear polarizations. Each 
brightness temperature is precisely calculable from thermcouple 
or thermistor measurements of Tbot *od Tcold » along with a- 
priori knowledge of the wire grid transmission and reflection co- 
efficients (6], the load emission and bistatic scattering properties, 
and the load background brightness temperatures. 

By rotating ths device about the feedhorn axis, the ra- 
diometer fesd can be illuminated by a set of known, linearly 
independent brightness temperature vectors. In practice, obser- 
vations would be made at four values of the load angle a: 0* , 
45 # , 90*, and 135*. The four [/-channel observations are related 
to the load temperatures by: 
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5. Expected Performance 

A linear estimator is used to infer 7 from the noisy mea- 


surements 


*/: 


f -*(#)• (*/-<*/» + <T) 


( 8 ) 


where {•) is the statistical mean. Using a minimum mean-squ are- 
error (MMSE) criterion, the optimum LMMSE estimator is: 

T * + ^*r»r] (®) 

where Bjt * the covariance matrix for the brightnses vector 
in the natural basis. The measurement error eo variance matrix 
becomes: 

3S„ * <(f-T)(*-T)‘> 

« (7-EF)*rr(7-EF)' + ( 10 ) 
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Fig. S shows LMMSE calculations of (he expected srrors 
in tv* fa* “d fv-a mfv -fa function of ^ for 4 typical 
conical scanner. The calculations assume #t * 1 K, » stan* 
dard deviation of 50 K for the vertical and borisontal brig fatness 
temperatures, and a standard deviation of 7 K for the vertical* 
faorisontaJ brightness difference Tv -a *Ty -Ta- 


1 


TT 


2500 2475.5 0 

2475.5 2500 0 

0 0 0 


( 11 ) 


These statistics are expected for conical swaths with — 45* half- 
cone angles [2]. The vertical and horisontal electric fields are 
assumed to be un correlated (Ter = 0), which is equivalent to the 
Miumption that the V — E basis is the natural basis. 


Calculations of the expected measurement error using the 
LMMSE estimator were performed for fixed-feed dual-polarisation 
radiometers configured both with ("5-parameter*) and without 
("2-parameter*) the -channel (Fig. 5). The electronic polarisa- 
tion basis rotation technique suggests significantly smaller mea- 
surement errors in fv, fa » end the difference fv-B than pos- 
sible using a fixed-feed dual-polarised radiometer without the 
(/•channel. The improvement is particularly striking for polar- 
isation skew angles near 45*. With the electronic technique, the 
measurement error is practically independent of 4* and hence, 
independent of the scan position within the conical swath. 


For a cross-track scanner, the antenna beam is oriented 
over a range of observation angles (typically from nadir+45* to 
nadir-45*) during the scan. Near nadir, the brightness difference 
Ty-jr drops to sero, (i.e., Tv and Tg become indistinguishable). 
The correlation between Ty and Tjr ae a function of observa- 
tion angle 9 (measured with respect to nadir) can be modelled 
approximately by: 


RTvT„ * 2475.5 + 24.5cos(20) (12) 


This relationship has been incorporated into the LMMSE estima- 
tor. Assuming that the nadiral scanner spot has a polarisation 
skew angle 4 « 45*, smaller measurement errors in fv~a result 
near 4 = 45* than in the case of tbe conical scanner (Fig. 4, 
"5-parameter*). However, since tbe nadiral Ty - Tg correla- 
tion can also be exploited using a conventional dual -polarised ra- 
diometer (with no (/-channel) and a LMMSE estimator (Fig. 4, 
"2-parameter*), a smaller reduction in error is realised by the 
addition of tbe 27-channel to a croee-track ecanner than when 
added to a conical scanner. 


6. Discussion 

Because of potential applications in conical and croee- 
track scanning microwave precipitation radiometers, we plan to 
prove the concept of electronic polarisation rotation through both 
laboratory and aircraft experiments. Our laboratory test will use 
the NASA Goddard Space Flight Center’s 22-GHz dual-polarised 
fixed-beam radiometer. Overall performance will be judged us- 
ing Brswster-angle observations of a water surface, which ex- 
hibits strong, predictable microwave polarisation characteristics. 
Pending favorable laboratory results, we have proposed to imple- 
ment the technique on the 85-GHs chsnnel of the NASA Mar- 
shall Space Flight Center’s Advanced Microwave Precipitation 
Radiometer (AMPR). The AMPR is a cross-track microwave 


scanner configured for operation on board tbe NASA ER-2 high- 
altitude aircraft with channels at 10, 12, 37, and 85 GH 2 . Al- 
though the error reduction attainable using electronic polariza- 
tion basis rotation on a cross-track scanner is not as large as that 
attainable on a conical scanner, it is nonetheless large enough to 
be useful. 
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Figure 2: Proposed dual-polarised calibration load for use with 
the 5- channel radiometer. The angle o is variable, thus produc- 
ing a set of known brightness vectors (T/A,T/tr,T/*). 
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Figure 1: Block diagram of a practical 3-channel (T/a,!/*, and 
T/v) radiometer to implement electronic polarization basis rota- 
tion. 


Conical Scanner: LMMSE Inversion X-Track ScannerrLMMSE Inversion 



Figure 3: Comparison of expected errors for a conical scanner em- 
ploying LMMSE polarization basis rotation without (2-parame- 
ter) and with (S- parameter) information from the ^-channel. 
Error estimate curves are shown for measurements of Ty (or 
equivalently, Tg) and Ty - Tg. 


Figure 4: Comparison of expected errors for s cross-track scanner 
employing LMMSE polarization basis rotation without (2-param- 
eter) and with (3-parameter) information from the 17-channel. 
Error estimates are shown for measurements of Ty (or equiva- 
lently, Tg) and Ty - 7*. The LMMSE estimator uses the corre- 
lation between TV and Tg at scan angles near nadir (i.e., skew 
angles near 45*). 
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APPENDIX B 


(Sufcroitted to the 1992 International Conference an Micrcwave 
Symposium, Albequerque, NM, June 1-5, 1992.) 

Laboratory Demonstration of 
Electronic Polarization Basis Rotation 
A.J. Gasiewski, D.B. Kunkee 
School of Electrical Engineering 
Georgia Institute of Technology 
Atlanta GA, 30332-0250 

Summary 

In microwave Earth remote sensing, it is often desirable to build 
passive imaging radiometers with dual orthogonal-linear polarization sen- 
sitivity. The additional information obtained by observing both vertically 
and horizontally polarized brightness temperatures can be used to im- 
prove estimates of rainfall, water vapor, and ocean surface winds, and 
to enhance the detectability of sea ice. However, in order to retain co- 
incidence between the radiometer feedhorn ’s polarization basis and the 
Earth’s vertical-horizontal basis, previous mechanically-scanned imager 
designs have required that the feedhorn and associated receiver electron- 
ics be physically rotated along with the scanning mirror. This is both 
cumbersome and expensive. 


By measuring the first three feedhorn Stokes parameters (Tx, Tg, 
and Tu), the measured Stokes vector can be transformed into any desired 
basis, including the natural basis defined by the local vertical and hori- 
zontal directions. The technique is called “Electronic Polarization Basis 
Rotation” [1], A practical implementation of EBPR was proposed by 
Gasiewski [2], whose scheme included a polarized calibration load. In this 
article, we report the first successful demonstration of a well-calibrated 
polarization correlation radiometer with EBPR capabilities. 


The radiometer (depicted schematically in Figure 1) is a 90-GHz 
dual orthogonal linear polarization unit with three outputs: one for each 
of the two orthogonal polarizations (t> a and vj,) and a cross-correlation 
channel (v u ). The v u channel is formed by an adding correlator with the 
addition of a post-detection summing circuit. This circuit removes the 
relatively large signals caused by the squaring of the IF voltages in each 
of the channels A and B. The result is three signals that are linearly 
related to the first three feedhorn Stokes’ parameters: 


+ V + IT (1) 

where V represents the three video output signals from the radiometer 
depicted in Figure 1, <7 and <5 are the instrument gain and offset pa- 


u = J 


T a 

Tu 

T b 
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rameters, and ?T is the instrument noise referred to the video outputs. 
Dicke-switching is used to reduce the effects of receiver gain and offset 
fluctuations. The phase difference of the local oscillator signal at each 
mixer was nulled using a phase shift A <p. In addition the path length 
of the two IF signals was equalized to well within a correlation length 
l e ~ Vp/W , where v p is the phase velocity, and W is the bandwidth 
(Figure 1). 

The parameters ^ and 7i are needed to convert the measured voltage 
data into feedhorn brightnesses. However, both ^ and 5 remain stable 
only over a time scale of minutes, and thus need to be periodically esti- 
mated. Moreover, interchannel crosstalk (indicated by nonzero diagonal 
elements in occurs. For example, imperfect adjustment of the sum- 
ming circuit causes the off-diagonal elements in the second row of g to 
be nonzero. Similarly, polarization crosstalk causes the off-diagonal ele- 
ments in the first and third rows of f to be nonzero. The measured cross 
polarization rejection was approximately 15 dB. 

A calibration procedure employing a unique polarized calibration 
load (Figure 2) was successfully used to estimate all elements in g and o. 
The procedure requires an observation of the calibration load assembly at 
three unique values of c* in addition to an observation of a non-polanzed 
load. Four different observations (each resulting in 3 measured values) are 
required to achieve an observation matrix of full rank enabling estimation 
of all twelve elements in g and o. 


Once p and o have been determined, the feedhorn brightness tem- 
peratures (T/) can be estimated from the video outputs, 



Ta 

Tv 

Tb. 




( 2 ) 


The value of 7 / is then used to estimate 7 , the three Stokes parameters 
in the Earth’s vertical-horizontal basis [1]: 
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( 3 ) 


where <f> is the polarization skew angle as shown in Figure 1. 


Near Brewster-angle observations of a water surface using a coni- 
cally scanned mirror demonstrate the ability to rotate the polarization 
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basis, that is to obtain 7 from Tj. A conical scan was performed with 
an observation angle of 65° with respect to vertical (nadir) at a, height 
of a few feet over the water surface. To continuously observe 7 during 
the scan, (3) was evaluated using values of 4> determined by the scan- 
ning mirror orientation with respect to the feedhorn basis. Under these 
conditions, 7 remains constant while 7/ changes significantly due to the 
polarizing nature of the water surface. The results of the experiment 
are shown in Figure 3. The predicted value of 7 based on radiative 
transfer calculations appears as the solid straight lines in Figure 3. A 
small systematic error (< 7 K) caused by mechanical misalignment of the 
scan mirror and feedhorn axis is noticeable in the graph of T v and 7*. 
Otherwise, rotation of the feedhorn brightnesses is straightforward; the 
rotated brightnesses track the calculated vertical and horizontal bright- 
nesses well. 
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Scan of Fresh Water Surface at 90 GHz 



Figure 3: Scan of a fresh water surface using the 3-channel radiometer. 
The parameters 7U, Tv and Tb are the measured feedhorn brightness 
temperatures; T v , 7* and U are estimates of the brightness temperatures 
in the Earth’s vertical and horizontal basis. The predicted values of T v 
and 7\ from a numerical radiative transfer model are illustrated by the 
heavy straight lines. 
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